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ABSTRACT 

We present Sunyaev-Zel'dovich (SZ) measurements of 15 massive X-ray selected galaxy clusters 
obtained with the South Pole Telescope (SPT). The SZ cluster signals are measured at 150 GHz, and 
concurrent 220 GHz data are used to reduce astrophysical contamination. Radial profiles are computed 
using a technique that takes into account the effects of the beams and filtering. In several clusters, sig- 
nificant SZ decrements are detected out to a substantial fraction of the virial radius. The profiles are 
fit to the j3- model and to a generalized Navarro- Frenk- White (NFW) pressure profile, and are scaled 
and stacked to probe their average behavior. We find model parameters that are consistent with pre- 
vious studies: f3 = 0.86 and r cor e/^500 = 0.20 for the /3-model, and (a n , /3 n , j n , C5oo)=(l-0, 5.5, 0.5, 1.0) 
for the generalized NFW model. Both models fit the SPT data comparably well, and both are con- 
sistent with the average SZ profile out to beyond rsoo- The integrated Compton-?/ parameter Y$z is 
computed for each cluster using both model-dependent and model-independent techniques, and the 
results are compared to X-ray estimates of cluster parameters. We find that Ysz scales with Yx and 
gas mass with low scatter. Since these observables have been found to scale with total mass, our 
results point to a tight mass-observable relation for the SPT cluster survey. 

Subject headings: cosmic background radiation - cosmology: observations - galaxies: clusters: general 
- galaxies: clusters: intracluster medium 
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1. INTRODUCTION 

Galaxy clusters are the largest known gravitationally 
collapsed objects. They are believed to have taken nearly 
a Hubble time to form, and their abundance as a function 
of redshift is crucially dependent on the composition and 
expansion history of the universe. A sufficiently large 
and well-understood sample of galaxy clusters can there- 
fore be used to constrain cosmological models. Several 
survey campaigns using a variety of techniques are cur- 
rently underway. Precision measurements of individual 
clusters can be used to study the mass-observable rela- 
tions of these surveys, as well as to study the structure 
of the intracluster medium (ICM). 

Well-established techniques for probing galaxy clusters 
include optical, radio, and infrared observations of mem- 
ber galaxies, gravitational lensing of background galax- 
ies, and X-ray measurements of bremsstrahlung from 
hot intra-cluster electron gas. The Sunyaev-Zel'dovich 
(SZ) effect, which occurs when cosmic microwave back- 
ground (CMB) photons inv erse Compton scatter off the 
hot electrons of the ICM (jSunvaev fc Zerdovichl 119721 : 
lBirkinshawlll999[). provides a com plementary and pow- 
erful probe (|Carlstroni 
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The SZ effect can be approximated as the sum of two 
components: that caused by the random thermal mo- 
tion of the hot scattering electrons, and that caused by 
the cluster peculiar velocity relative to the Hubble flow. 
The former, known as the thermal SZ effect, leads to a 
distortion of the CMB Planck spectrum of the form 
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4)(l + 5sz(x,T e )), (1) 



where x = hv /ksTcMB is the dimensionless frequency, 
Tcmb is the temperature of the CMB, T e is the tem- 
perature of the ICM electrons, and Ssz is a relativistic 
correction 22 . For nonrelativistic electrons, the distortion 
results in a decrement in the measured CMB temperature 
at frequencies below ~ 217 GHz and an increment above. 
The amplitude of the effect, expressed as a change in tem- 
perature AT relative to the CMB temperature T CMB , is 
given by 
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where y is the Compton-y parameter, a measure of inte- 
grated pressure; gt is the Thomson cross section; n e and 
T e are the electron number density and temperature; and 
the integral is along the line of sight through the clus- 
ter. The integrated y parameter out to a given radius 
is a quantity proportional to the total thermal energy 
of the cluster, and is expected to provide a low scatter 
estimate of the total cluster mass ([da Silva 20(3; iNagail 
120061 : iMotl et al.ll2005t iKravtsov et al.ll2006l ) . For typical 
massive clusters, the kinetic SZ effect produces a signal 
that is much smaller than the thermal SZ signal except 
at frequencies near the ~ 217 GHz thermal SZ null. 

Mu ch of the past work on the thermal SZ effect 
(e.g.. iGrego et all 2000; Reese et al.l 120021 : iBenson et all 
I2QQ4J : lHalverson et al.l l2QQ9h has modeled cluster SZ sig- 
nals by the 0— model ( Cavalier e fc Fusco-Femiano1ll976l : 
I Cavalier e fc Fusco-Femian o 1978). The functional form 
of this model can be derived from a parameteriza- 
tion of density under the assumption of isothermal- 
ity, but since isothermality is a poor assumption for 
many clusters — particularly for merging and cooling core 
systems — we instead consider the f3— model simply as 
a fitting function without a strong physical motiva- 
tion. Measurements and simulations show significant 
departures from this model within the cluster core and 
outside ~ r25oo 7 the radius at which the cluster mat- 
ter density d rops to 2500 times the critical density of 
the universe ( Vikhlini n et al.ll2QQ5t iPiffaretti et al.ll2QQ5t 
lHallman et al.ll2007[ ). lLaRoque et al.l (|2006l ) examined 
two extensions of the j3— model using a sample of 38 
clusters with both X-ray and SZ data. As an alter- 
native to the /?— model, Nagai e t al.l (l2007j here after 
NO 7) have proposed , and Mr oczkowski et al.l ([2009) and 
lArnaud et al.l (j2QQ9h have investigated, a generalization 
of the Navarro- Frenk- White (NFW) model as a parame- 
terization of pressure in clusters. This model is expected 
to describe the cluster pressure out to a significant frac- 
tion of the virial radius. 

In this work, we describe a set of deep SZ cluster obser- 
vations at 150 and 220 GHz undertaken with the South 
Pole Telescope (SPT) during the austral winter of 2008. 

22 We adopt the relativistic correction from Itoh et al. ( 2000). 



Data from the 150 GHz band are used to measure the 
thermal SZ signal, and concurrent data from the 220 GHz 
band — which is near the thermal SZ null — are used to 
reduce the effect of astrophysical contamination. The 
220 GHz data are too noisy to allow a detection of the 
kinetic SZ effect, so in this paper we focus exclusively on 
the thermal SZ effect. 

Of the 13 highes t-luminosity REFLEX clusters 
(Boh ringer et al.l l2004f ) in the range of elevation angles 
accessible to the SPT, 11 make up the primary sample 
for this work; a supplemental sample of four additional 
clusters is also included. Both samples are listed in Ta- 
ble [U Taking advantage of the SPT's wide field imaging 
capabilities, we make a model independent estimate of 
the radial profiles of each cluster and place constraints 
on the integrated pressure out to large radius. We also 
fit the profiles both to the (3 — mod el and to the general- 
ized NFW model proposed bv lN07l (hereafter the GNFW 
model). Using a stacked analysis, we estimate the slope 
of the pressure profile in the cluster outskirts. Finally, 
we estimate the integrated y parameters for each cluster 
and compare to X-ray results. 

2. INSTRUMENT, OBSERVATIONS, AND DATA 
REDUCTION 

2.1. Instrument 

The SPT is an off-axis Gregorian telescope with a 10 
meter diameter primary dish. O ne of a new gener ation 
of SZ instr uments (such as ACT ([Hincks et alJ l2QQ9) and 
APEX-SZ (jHalverson et aLll2QQ9[ |). the SPT has several 
advantages for cluster observations: 

• Its resolution (~ V at 150 GHz) is matched to clus- 
ter scales. 

• It is able to map clusters efficiently due to its high 
sensitivity over a wide field of view (~ 1 deg 2 ). 

• It observes concurrently in multiple bands, includ- 
ing one band centered near the 217 GHz SZ null, 
allowing primary CMB anisotropics to be separated 
from SZ signals. 

• The atmosphere at its observing site, roughly 1 km 
from the geographic South Pole , is the best on 
Earth for millimeter astronomy (Bussm ann et al.l 
l2005h . 

The SPT is optimized to image large areas of the CMB 
sky to arcminute resolution with a primary goal of 
identifying a sample of massive clusters out to high 
redshift. The fi r st su ch clusters were reported in 
IStaniszewski et al.l (|2009L hereafter S09). 

The SPT receiver is equipped with a 960-element array 
of superconducting transition edge sensor (TES) bolome- 
ters, each read out by a frequency-domain- multiplexed 
system using superconducting quantum interference de- 
vice (SQUID) amplifiers. This detector array is divided 
into six triangular "wedges," each of which is populated 
with detectors sensitive to radiation within a single fre- 
quency band. During the 2008 observing season, the 
band centers were ~95 GHz for one wedge, ~150 GHz 
for three wedges, and ~220 GHz for the remaining two. 
Due to a defect in detector fabrication, few of the 95 GHz 
detectors were active during the 2008 observing sea- 
son, and so only the two higher-frequency bands are 
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included in this work. The bands were measured with 
a Fourier transform spectrometer, and the bands were 
found to be Av = 35.6 GHz at v — 152.6 GHz and 
Av = 42.7 GHz at v = 220.0 GHz. In a typical obser- 
vation used in this paper, the instantaneous per-detector 
sensitivity was ~ 500/iK CMBV /s for the 150 GHz detec- 
tors, while the 220 GHz detectors had sensitivities rang- 
ing from - 1000 - 1500/iK nMRA /s 23 . For further de - 
tails on the telesco pe an d receiver. see | Ruhl et all ([2004); 
iPadin et al.l (|2008f ). and ICarlstrom et al.l (|2009f ). 

2.2. Observations 

The clusters discussed in this work were observed us- 
ing constant-elevation scans, which entail sweeping the 
telescope at constant angular velocity in azimuth across 
the field and back, stepping in elevation, and repeat- 
in g. Th is observing strategy is similar to that described 
in IS09L but over a much smaller field — approximately 1 
deg 2 . The azimuthal scan speed is set such that ~ 0°.25 
on the sky is scanned in 1 s, and each elevation step is 
0°.005. A complete set of scans across the entire region, 
which we refer to as one observation, takes approximately 
55 minutes to complete. The observations are then com- 
bined to produce a single map in each band for each field. 
The inner region of each 150 GHz map, approximately 
0°.5 x 0°.5 on the sky, has even coverage to within 10%. 
Coverage is less uniform at 220 GHz due to the spatial 
distribution of the detectors on the focal plane. 

One of the clusters in the primary sam ple, A S 0520, 
overlaps with the survey field described in IS09I . We re- 
processed the 2008 150 GHz and 220 GHz survey data 
with filtering equivalent to that described in Section [2^ 
and split the survey data into subsets, each with a weight 
per pixel equal to that of the targeted observations. This 
allows the survey data to be treated in a similar manner 
to the data from the targeted observations for the pur- 
poses of calculating noise estimates. Adding the survey 
data effectively doubles the number of targeted observa- 
tions of this cluster. 

Between individual observations, we perform a series 
of shor t calibration measurem ents described in more de- 
tail in ICarlstrom et aTl ([2009) . These include measure- 
ments of a chopped thermal source, ~ 2° elevation nods, 
and scans across the galactic HII regions RCW38 and 
MAT5a. This series of regular measurements allows us 
to identify detectors with good performance, assess rela- 
tive detector gains, and monitor atmospheric opacity. 

2.3. Data reduction 

The data reduction process used for this work is similar 
to the preliminary survey pipeline described inlS09], with 
some minor differences noted below. 

2.3.1. Data selection and filtering 

The first step in the processing pipeline is to identify 
the data that will be included in each single-observation 
map. For every observation, a set of well-performing 
detectors is identifi ed b y applying a set of noise and 
sensitivity criteria (|S09l ). On average, 302 detectors 
at 150 GHz and 165 at 220 GHz pass this initial cut. 

23 Throughout this work, K C mb refers to equivalent CMB fluc- 
tuation temperature, i.e., intensity divided by dB u /dT (at T = 
2.73K). 



The data are then divided into individual azimuth scans. 
0.7% of scans have large instantaneous pointing offsets 
or data acquisition problems, and are omitted from fur- 
ther processing. A further set of selection criteria is then 
applied to detectors on a scan-by-scan basis (excluding, 
for example, data contaminated by cosmic rays). Data 
from detectors passing these cuts, an average of 291 de- 
tectors at 150 GHz and 162 at 220 GHz for each scan, 
are included in the final maps. 

The receiver exhibits some sensitivity to the pulse- 
tube cooler used to cool it to ~ 4 K, resulting in oc- 
casional lines in the detector noise power spectra, so a 
small amount of bandwidth (< 0.4%) is notch filtered 
out of the detector timestreams. The timestreams are 
further processed by deconvolving the effect of the de- 
tector time constants. We also low-pass filter at 25 Hz, 
well above the signal bandwidth. 

For every scan, each detector's timestream is then fit 
simultaneously to a number of template functions, and 
the best fit to each template is subtracted out. This time- 
domain filter reduces the effect of low-frequency noise in 
the detector timestreams due to readout noise or to at- 
mospheric fluctuations. For the analysis presented here, 
we use as our templates a second-order Legendre polyno- 
mial and a common mode constructed from the wedge- 
averaged timestreams, resulting in a characteristic filter 
scale of roughly one half degree. This light filtering is suf- 
ficient to remove much of the low- frequency noise, due in 
part to the outstanding atmospheric conditions at the 
South Pole. 

The common mode template is constructed from the 
mean across all of the well-performing detectors in a 
given wedge using a nominal relative calibration. In IS09I . 
by contrast, the common mode template was constructed 
from the mean across the entire array for each band. 
Since the 2008 SPT receiver contains three 150 GHz and 
two 220 GHz wedges, the detectors in the two bands have 
significantly different distributions on the focal plane; 
thus the array common mode removes signals on differ- 
ent spatial scales for each band. Since the common mode 
templates for individual wedges are relatively similar to 
one another, the wedge common mode preserves signals 
on similar spatial scales for each band, simplifying com- 
parisons between the 150 GHz and 220 GHz maps. 

2.3.2. Calibration and beams 

We take advantage of the superb WMAP5 absolute 
calibration and calibrate the SPT data by comparing 
the 15 GHz SPT maps wi th WMAP5 V- and W-band 
maps (jHinshaw et al.ll2009f ). The SPT performed dedi- 
cated calibration scans of five large fields totaling 1250 
deg 2 of sky during 2008. The WMAP5 maps were then 
resampled with the SPT pointing information, and the 
resulting timestreams were passed through a simulated 
SPT analysis pipeline to capture the filtering applied 
during the SPT map-making algorithm. The ratio of the 
cross-spectra of the SPT maps and the filtered WMAP 
maps after correcting for the instrumental beams pro- 
vides an estimate of the relative calibration factor be- 
tween the two experiments. This procedure, similar to 
that used for the analy sis of ACBAR data and detailed in 
iReichardt et al.l (|2009l ) , results in an absolute calibration 
accurate to 3.6% at 150 GHz. 

The 150 GHz calibration is transferred to 220 GHz by 
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TABLE 1 

Summary of cluster sample. 



ID 




R.A. a 




Decl. a 








REFLEX 




Mgas 


Ref. fa 




















(10 44 erg s- 1 ) 


(keV) 


(10 








Primary sample 


A 2744 




h 14 m 18 s 


6 


-30°23 


15" 


4 


0.307 


12.92 


10.1 ± 0.3 


1.0 


± 0.2 


Zhang et al. 


(2006) 


RXCJ0217.2 


5244 


2 h 17 m 12 s 


6 


-52°44 


49" 


2 


0.343 


12.03 


10. 9 C 






Bohringer et 


(20(M ) 


RXCJ0232.2 


4420 


2 h 32 m 18 s 


8 


-44°20 


5l" 




0.284 


9.65 


7.0 ± 0.3 


0.9 


± 0.2 


Zhang et al. 


(2006) 


AS 0520 




5 h 16 m 35 s 


2 


-54°30 


36" 


8 


0.294 


13.87 


7.5 ± 0.3 


0.8 


± 0.2 


Zhang et al. 


(2006) 


RXCJ0528.9 


3927 


5 h 28 m 52 s 


5 


-39°28 


16" 


7 


0.284 


13.12 


7.2 ± 0.4 


0.9 


± 0.1 


Zhan^^L 


(2006) 


AS 0592 




6 h 38 m 46 s 


5 


-53°58 


18" 





0.222 


10.62 


8.0 ± 0.4 






Hughe^^l 


(2009 ) 


A 3404 




6 h 45 m 30 s 





-54° 13 


42" 


1 


0.164 


7.36 


8.1 ± 0.3 


0.9 


± 0.1 


Zhang et al. 


(2008) 


1ES 0657-56 
RXCJ2031.8 


4037 


6 h 58 m 30 s 
20 h 31 m 51 


2 
s .5 


-55°56 
-40°37 


33" 
14" 


7 



0.297 
0.342 


23.03 
12.04 


10.6 ± 0.2 
10.9° 


1.8 


± 0.3 


Z h a n ge^^ 
Bohringer et 


(2006) 
al. f2004l 


A 3888 




22 h 34 m 27 


s .1 


-37°44 


7". 5 


0.151 


7.31 


7.8 ± 0.4 


0.8 


± 0.1 


Zhang et al. 


(2008) 


AS 1063 




22 h 48 m 44 




-44°31 


44" 


4 


0.346 


30.79 


11.1 ± 1.1 


1.2 


± 0.1 


Maughan et 


al. (20081 








Supple 


mental sample 








RXCJ0336.3 


4037 


3 h 36 m 18 s 


7 


-40°37 


20" 





0.173 


2.53 


7.2° 






Bohringer et 


al. C20041 


RXCJ0532.9 


3701 


5 h 32 m 55 s 


9 


-37°l / 34 // .5 


0.275 


6.94 


9.5 ± 0.4 


0.6 


± 0.1 


Zhang et al. 


(2006) 


MACSJ0553. 4-3342 


5 h 53 m 27 s 


2 


-33°42 


53" 





0.407 




13 1 + 3 - 8 






Cavagnolo et 


al. C20081 


A 3856 




22 h 18 m 39 


S g 


-38°53 


43" 


6 


0.141 


3.78 


6.7 ± 0.2 


0.5 


± 0.1 


Zhang et al. 


(20081 



a R.A. and decl. are the right ascension and declination angles of the cluster X-ray centroid. 

b R.A., decl., redshift, and electron temperature T e are taken from these papers and references therein. REFLEX luminosities are from 

Bohr inger etaD (|200l ). 

c For the clusters without published X-ray temperatures, T e was estimated using a scaling relation from luminosity ( Markevitch et al. 1998). 
d Luminosity in rest frame 0.1-2.4 keV band corrected for any missi ng flux. 

e Temperatures are av erages over a range of radii: 0.1-0.5r5oo in Zhang et al. (2006), 0.2-0.5rsoo in Zhang et al. (2008 ), r < rspp in 
Maughan et al. (2008), and r < r25oo with core excised Cavagnolo et al. (2008). The radial range is unspecified in Hughes et al. (2009). 



TABLE 2 

Summary of SPT cluster maps. 



ID 


R.A. a 


Decl. a 


Observing Time 


150 GHz Depth b 


220 GHz Depth 6 








(hr) 


OK CMR ) 


OK^mr) 


Primary sample 


A 2744 


h 14 m 16 s .8 


-30 o 23'22".2 


6.83 


12.3 


36.2 


RXCJ0217. 2-5244 


2 h 17 m 13 s .3 


-52°44 / 49 // .5 


8.32 


11.5 


40.8 


RXCJ0232. 2-4420 


2 h 32 m 16 s .5 


-44°21 / 15 // .5 


10.46 


9.2 


32.6 


AS 0520 


5 h 16 m 35 s .4 


-54°30 / 19 // .3 


18.48° 


10.3 


17.4 


RXCJ0528. 9-3927 


5 h 28 m 52 s .0 


-39°27 / 57 // .6 


5.23 


14.1 


39.4 


AS 0592 


6 h 38 m 46 s .l 


-53°58 / 3l".3 


8.32 


18.9 


70.7 


A 3404 


6 h 45 m 29 s .3 


-54°13 / 20".0 


5.54 


14.5 


46.5 


1ES 0657-56 


6 h 58 m 29 s .7 


-55°56 / 38".2 


12.59 




22.5 


RXCJ2031. 8-4037 


20 h 31 m 51 s .l 


-40°37 / 23 // .2 


5.23 


16.3 


49.7 


A 3888 


22 h 34 m 26 s .6 


-37°44 / 57 // .9 


5.23 


15.4 


42.6 


AS 1063 


22 h 48 m 44 s .9 


-44°31 / 43 // .l 


12.21 


10.3 


29.7 


Supplemental sample 


RXCJ0336. 3-4037 


3 h 36 m 16 s 5 


-40°37'26".6 


5.23 


13.1 


48.1 


RXCJ0532. 9-3701 


5 h 32 m 55 s .l 


-37°l / 40".8 


10.46 


12.8 


38.1 


MACSJ0553. 4-3342 


5 h 53 m 24 s .5 


-33°42 / 35 // .4 


5.82 


15.4 


51.0 


A 3856 


22 h 18 m 37 s 6 


-38°53 / 54 // .8 


8.72 


11.0 


34.9 



a R.A. and decl. are the right ascension and declination angles where the average SZ decrement within a l f radius is the largest. 
b Depth is determined by the standard deviation of the jackknife maps smoothed with a Gaussian with FWHM=1 ; . 

C AS 0520 maps include data from the 2008 5 h 30 m survey, which effectively doubles the amount of targeted observing time. This value 
includes the factor of ~ 2. 



making use of the overlapping coverage on SPT's high 
signal to noise maps of CMB fluctuations in large sur- 
vey fields. We calculate the ratio of the cross-spectra of 
the map at each frequency band to the auto-spectra of 
the CMB-dominated 150 GHz map after correcting for 
the beam and filtering differences. The calibration uncer- 
tainty of the 220 GHz maps is found to be 7.2%. Because 
the thermal SZ signal in the 220 GHz map is small, the 
higher calibration uncertainty only marginally increases 
the errors on the cluster profiles discussed in this paper. 

The 2008 SPT beams were measured by observing 
bright sources and averaging the response for all detec- 
tors in a given band. The sidelobes of the beam were es- 
timated from observations of planets, and — since the de- 
tectors show signs of nonlinearity when pointed directly 
at planets — the main lobes were estimated from obser- 



vations of quasars. The main lobes are reasonably ap- 
proximated by two-dimensional elliptical Gaussians with 
average full widths at half maximum (FWHM) of 1.2 and 
1.1 arcmin for 150 and 220 GHz. These are larger than 
would be naively calculated from the aperture diameter 
and observing wavelength due to the combination of an 
underilluminated primary (to red uce spillover ) and tele- 
scope pointing uncertainty; see iPadin et al.l (j2QQ8h for 
details. The sidelobes fall off from ~ — 30dB at a radius 
of 5' to < — 45dB at 30', and account for about 15% of 
the degree-scale beam response. Errors in this beam es- 
timate contribute to the overall calibration uncertainty, 
and also add small (~ 1%) additional uncertainties to 
the best-fit model parameters. 

We model the response of the instrument to an as- 
tronomical source by constructing a simulated map of 
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a point source, convolving it with the beam, construct- 
ing simulated timestreams using telescope pointing infor- 
mation, and then running the timestreams through our 
processing pipeline. The output map normalized to the 
input signal is referred to as the transfer function, and 
is used to subtract bright point sources from the clus- 
ter maps. Representative transfer functions are shown 
in Figure [TJ The negative stripe at constant elevation 
is due to the polynomial template subtraction, and re- 
sults in apparent temperature increments to the east and 
west of the cluster proportional in amplitude to the SZ 
temperature decrement. This effect can be ameliorated 
by masking the cluster before applying the time-domain 
filtering (see Figure [2]), but doing so reduces the effective- 
ness of the atmospheric noise removal. Since we account 
for the effects of the time-domain filtering on the profiles 
as described in Section \4~T\ we avoid degrading the signal 
to noise of our results and analyze the unmasked data. 

2.3.3. Mapmaking and astrophysical background cleaning 

For every observation, a map is made for each ob- 
serving frequency using the processed data for all de- 
tectors in that band. Pointing information (R.A. and 
decl.) is calculated for each detector using focal-plane 
offsets measured in observations of the galactic HII re- 
gions, and boresight pointing calculated using data from 
the telescope pointing readout system, with a set of 
corrections described in S09 that result in a point- 
ing rms of 7" . These coordinates are then converted 
to pixel number using a S anson-Flamsteed projection 
(|Calabretta fc Gre isen 20Q2j) with 0.25 arcmin pixels. All 
measurements of a given pixel's brightness are averaged 
using inverse-variance weighting based on the mean of 
each detector's processed and relative-gain-scaled power 
spectrum between 1 and 3 Hz, which for our scan speed is 
a frequency range well matched to the scale of SZ cluster 
signals. 

All of the single-observation maps for a given cluster 
in a given band are then co-added to produce a set of 
single-band maps, which we denote as M 150 and M 220 
for the 150 GHz and 220 GHz bands. Each single- 
band map consists of signals from several astrophysical 
sources: the SZ effect, which is at arcminute scales and is 
much stronger at 150 GHz; the primary CMB anisotropy, 
which dominates at large spatial scales and is the same 
in both bands; and a background of typically faint dusty 
point sources, which are at small angular scales and are 
stronger at 220 GHz. We wish to produce a "band sub- 
tracted map" M sub , in which the 220 GHz map is used 
to remove a fraction of the CMB and point source back- 
ground while leaving the SZ signal intact. Since the 2008 
SPT receiver was less sensitive to CMB fluctuations at 
220 GHz than at 150 GHz, and since the atmospheric 
noise is worse in the higher-frequency band, this back- 
ground removal must be handled carefully to avoid in- 
troducing excess noise into the band subtracted map. 

Using our knowledge of the spatial scales of the CMB 
and point source signals, we apply a matched spatial fil- 
ter ip to the 220 GHz map, which we construct by re- 
quiring that t he variance in the band subtracted map 
be min i mized (Haehnelt & Tegmark 1996; Herranz et al.l 
[200230; iMelin et al.ll2006h . Since the signal and noise 
are more easily characterized in the spatial frequency do- 
main, we adopt the flat sky approximation and construct 



the filter as a function of multipole moment £. Denot- 
ing the Fourier transform of X as X, we can express the 
value of the band subtracted map at a given multipole 
moment as 

Mp b = M| 50 - ^Mj 20 . (3) 
The variance at a given t, T/g, is then given by 
V t = {\ -fa) 

( S? MB + SP + Nj 50 S? MB + aSf s \ 
V Sf MB + aSP S™ B + a 2 SP + N? ) 

(-k)' 

where S^ MB is the CMB signal, S^ s is the point source 
signal, is the noise in band X, and a is a factor 
corresponding to the spectral index of the point sources. 
The filter that minimizes this variance is given by 

gCMB + a gPS 

^ S™ B + a2Sp + Np°' [ ' 

Note that both the signal and the noise terms are 
anisotropic with respect to azimuth and elevation due 
to the scan strategy employed in these observations, ne- 
cessitating a two-dimensional filter function. Errors in 
either the signal or noise terms will result in a slightly 
sub-optimal filter, and will thus increase the errors in the 
results, but will not lead to systematic misestimations of 
the profiles or cluster parameters. 

The signal covariances of the primary CMB 
anisotropics and undetected point sources are esti- 
mated using simulations b ased on the best-f it WMAP5 
CMB powe r spe ctrum ([Nolta et al.l l2QQ9h and the 
iBorvs et"aTI (|2003l ) model for dusty point sources. We 
assume that the point sources are Poisson distributed 
and have a spectral index of 2.7 24 . The power spectrum 
of the source co unt distribution is c omput ed using 
the formalism in I White fc Majumdarl (|2004), and is 
inflated by a factor of 40% to account for le nsing by 
our massive cluster targets ([Lima et al.l l2009f ). Even 
with this enhancement, the dusty sources are predicted 
to contribute an rms of ~4 fiK to the band subtracted 
maps, well below the map noise level. We disregard 
contamination due to radio sources, as explicit simula- 
tions demonstrate ([Sehgal et al.ll2010l ) that they are not 
likely to significantly fill in the 150 GHz SZ decrements 
of massive z ^ 0.4 clusters, and as in most cases we are 
able to remove bright sources from existing catalogs (see 
Section 13. ip . For both the CMB and the dusty point 
sources, we generate 300 simulated maps and apply 
the same time-domain filtering that was applied to the 
SPT maps. We find 6^ MB and S% s by taking the mean 
of the two-dimensional spatial power spectra over the 
simulated maps. 

The instrumental and atmospheric noise properties 
of the maps are estimated using the two-dime n sional 
power spectra of noise maps ([Savers et al.l [2009; 

24 SPT data indicate a steeper spectral index between 150 and 
220 GHz (jHall et al.|[2009h , but we adopt this value— extrapolated 
from higher-frequency measurements — so that our estimate of 
point source power is an upper limit. Our results are insensitive to 
this choice. 
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Fig. 1. — Point source transfer functions at 150 GHz (left) and 220 GHz (right). Units are percent of maximum. The stripe through the 
center of the maps is the result of the polynomial template removal. 
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Fig. 2. — Maps of 1ES 0657-56 at 150 GHz, with the cluster region unmasked (left) and masked (right) for time-domain filtering. Units 
are /iXcMB • I n the masked map, data from a circular region 8' in radius centered on the cluster are not included in the polynomial and 
wedge common mode construction. This map gives a clear picture of the shape and amplitude of the SZ emission, but has slightly higher 
noise. The artifacts of the time-domain filtering in the unmasked map are taken into account in the cluster profiles by the algorithm 
described in Section [4~T1 



lHalverson et al. 2009). Under the assumption of station- 
arity in the map basis, the noise covariance matrix in 
band X, iV^ ise , is diagonal in the spatial frequency do- 
main and equal to the noise power spectrum. For each 
observation, we produce one map using only data from 
the left-going scans, and one map using only data from 
the right-going scans. We then multiply one half of the 
2n maps by —1 and co-add to produce a jackknife noise 
map. We repeat this process for m » n combinations 
of observations, computing the two-dimensional spatial 
power spectrum for each individual jackknife map. The 
a verage of thes e powe r spec tra is our estima t e of N * oise . 

Sav ers et al.l (|2009f ) and lHalverson et al.l (|2009f ) use 
this noise estimate to construct a covariance matrix and 
to fit models directly to their cluster maps. However, 



since we have relatively few independent observations of 
each cluster, our noise estimate is less well-constrained. 
We cannot combine data from different clusters to im- 
prove the noise estimate due to the fact that the noise is 
non- stationary between the observations. The amplitude 
of the atmospheric noise varies with time and as a func- 
tion of elevation angle, and the clusters in this sample 
are at a wide range of elevation angles. We circumvent 
this limitation by computing the projected radial profiles 
of the clusters, for which the covariance matrix has fewer 
degrees of freedom and is thus better constrained by our 
data. 

3. SAMPLE SELECTION AND CLUSTER MAPS 

During the austral winter of 2008, the SPT observed 
11 of the 13 highest-luminosity clusters in the REFLEX 
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Fig. 3. 150 GHz 1ES 0657-56 masked map with X-ray plasma 
density (white) and weak lensing surface mass density (black) con- 
tours overlaid. Units are fiKcMB- This cluster is an ongoing 
merger, and the electron gas is not in equilibrium within the grav- 
itational potential well. The SZ signal tracks the X-ray plasma 
density more closely than the surface mass density. 



survey within the range of declination angles observ- 
able by the SPT (-30° > decl. > - 70°). These clusters 
form the primary sample used in this work. One of the 
two REFLEX clusters that were omitted, RXCJ1253.6- 
3931, is associated with a radio bright PMN source that 
would have been difficult to remove from the SZ maps. 
The other, RXCJ1234. 2-3856, was omitted due to time 
constraints. Of the next three brightest REFLEX clus- 
ters, RXCJ0532. 9-3927 is a part of t he supplemental 
sampl e; RXCJ201 1.3-4037 was found by lMaughan et al.l 
(2008) to have an electron temperature of just 3.8 keV, 
and was not detected at greater than 3<r by the SPT; 
and RXCJ1317. 1-3821 is associated with a bright radio 
source. From the remaining 2008 SPT targeted clus- 
ter observations, four additional clusters were selected to 
form our supplemental sample. These clusters were cho- 
sen because of the availability of archived X-ray satellite 
data, which will be relevant for future work. They also 
contribute to our understanding of pressure profiles, and 
demonstrate the ability of the SPT to make high signal 
to noise maps of clusters with lower X-ray luminosities. 

The positions and depths of the final co-added maps 
are listed in Table [2j and the images can be found in the 
Appendix. Four maps are shown for each cluster: the 
150 and 220 GHz single-band maps, the band subtracted 
map, and a jackknife band subtracted map computed 
by multiplying half of the observations by —1 before co- 
adding. Each map is smoothed by a Gaussian with an 
FWHM of 1'. All clusters are detected with high signif- 
icance. Since the maps each have distinct noise proper- 
ties, we do not expect the detection significance to pro- 
vide quantitative information about cluster parameters. 

As discussed in Section I2.3.2| the time-domain pro- 
cessing creates large-scale distortions in the maps visible 
as temperature increments to the east and west of the 
clusters. No mask has been applied to the data used to 
construct the maps in the Appendix. Figure [3] shows the 
1ES 0657-56 SZ map constructed from data with the clus- 
ter masked before time-domain filtering. Shown as con- 
tours in Figure [3] are publicly available X-ray plasma sur- 



face density and weak lensing surface density data 25 . As 
expected, once the effects of the beam and time-domain 
filtering are mitigated, the SZ signal closely tracks the 
electron plasma density. 

3.1. Source contamination 

We attempt to remove bright point sources from the 
SPT maps. For sources more than 2' away from the clus- 
ter center, we identify and remove all sources detected at 
5cr or greater in either frequency band. A point source 
within 10" of quasar J001341. 2-300926 was detected at 
> 5<j significance in the map for A 2744, 16'. 5 from the 
cluster center, and has been removed from the maps by 
subtracting the point source transfer function scaled to 
the best-fit amplitude. The SPT measured fluxes for 
this source are 17.7 mJy at 150 GHz and 27.8 mJy at 
220 GHz. In the AS 0592 map, a point source corre- 
sponding to the radio source SUMSS J063845-535824 was 
removed. This source is ~ 1' from the cluster center, and 
the fluxes at 150 GHz and 220 GHz are 21.8 mJy and 
13.4 mJy, respectively. 

Sources within 2' of the cluster center may fill in 
the cluster SZ decrement without being directly de- 
tected. We searched existing radio source catalogs 
within 2' of our cluster locations using the NED 
database 26 . Two SUMSS sources with 843 MHz fluxes 
greater than 10 mJy were identified: SUMSS J021714- 
524529, a 24.5 mJy source near RXCJ0217.2-5244; 
and SUMSS J0 65837-555718, a 79.4 mJy source near 
1ES 0657-56. IColafrancesco et al.l ([2007ft note that 
the latter source has been found to have a power law 
SED with a = —0.9, and would therefore have fluxes 
at 150 GHz and 220 GHz that are well below our 
noise levels. Insufficient data exist to extrapolate the 
SUMSS J021714-524529 flux to 150 GHz. We assume 
that its flux contributes negligibly to the SZ flux at 150 
GHz, as would be the case for typical radio source SEDs. 
We note that RXCJ0217.2-5244 is not included in our 
scaling relation analysis. 

Dusty galaxies m ay also appear as brig ht point sources 
in the SPT maps. IWilson et all (|2008h report a dusty 
point source approximately V to the east of 1ES 0657-56, 
with flux density of 13.5 ± 0.5 mJy at 270 GHz. Assum- 
ing a spectral index of 2.7, we would expect to see 7.8 
and 2.8 mJy in the 220 and 150 GHz SPT maps. The 
SPT does not make a significant detection of a source 
at these coordinates in either band. We measure a tem- 
perature increment of 34 ± 33 /xKcmb in our 220 GHz 
map (smoothed by a V FWHM Gaussian). By finding a 
best-fit amplitude for our point source transfer function, 
we infer a 220 GHz flux density of 4.2 ± 4.5 mJy. At 
150 GHz, the strong SZ decrement near the source po- 
sition prevents us from determining the source flux. We 
find from simulations that a source with the predicted 
fluxes at these coordinates would have a < 5% effect on 
the corresponding bins in the band subtracted cluster 
profile. 

4. PROJECTED RADIAL PROFILES 
4.1. Profile computation 

25 ht tp://flamingos.astro.ufl.edu/le0657/public.html 

(|Clowe et al.l[2006T ) 

26 http://nedwww.ipac.caltech.edu/ 
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To proceed with measurements of global cluster prop- 
erties from the maps, it is useful to approximate the clus- 
ters as azimuthally symmetric and to calculate their pro- 
jected radial profiles, which are expected to take a uni- 
versal form ([Arnaud et al .1120091 ) . We do so by computing 
transfer functions for a set of concentric annuli with in- 
ner and outer radii set to the boundaries of the radial 
bins for which we wish to estimate an average SZ decre- 
ment. These transfer functions are then used to compute 
a profile transfer function matrix, which is inverted and 
applied to the maps in order to account for the effects 
of the beam and the time-domain filtering. This tech- 
nique also allows us to characterize the correlated errors 
between bins. 

The cluster center coordinates must be determined in 
order to generate radial profiles. We define the cen- 
ter position by fitting a /3-model with j3 = 0.86 to the 
map with the center R.A. and decl., central tempera- 
ture decrement, and core radius free to vary. We then 
marginalize over the central decrement and core radius 
to determine the cluster center. The fit uncertainties 
on the center position vary from 1" to 4", indicating 
very high signal to noise detections, and the results are 
consistent when the generalized NFW model is used in 
place of the j3— model. The difference between the X-ray 
centroid (which depends primarily on density) and the 
SZ center position (which depends on integrated pres- 
sure) is small for most clusters, but reaches 40 — 50" 
for RXCJ0232.2-4420, MACSJ0553.4-3342, and A 3888. 
An X-ra y-derived pressure map f or RXCJ0232. 2-4420 is 
shown in lFinoguenov et al.l ([2005f h and the SPT centroid 
falls within the elongated central high-pressure region. 

For each band subtracted cluster map M sub , we divide 
the 0°.5 x 0°.5 region surrounding the center position into 
radial bins b{. The bin spacing is defined so as to measure 
the profile of the cluster center with beam-scale resolu- 
tion, and to keep the number of bins small enough so 
that the noise properties can be well constrained. Each 
bin is then assigned a radius T{ by finding the average 
radius over all pixels within the bin. For bin bi and 
band X, we define Ti to be the mean SZ temperature 
decrement on the sky within the bin. We construct a 
set of simulated maps corresponding to each bin in each 
band, in which the "temperature" for map i is set to 1 
inside bin bi and elsewhere. These maps are then con- 
volved with the appropriate beam. To account for the 
effect of the time-domain filtering, we construct simu- 
lated timestreams from these maps using the telescope 
pointing data from each observation, and apply to them 
the same polynomial and common mode filtering that 
was applied to the data. We then produce from these 
filtered timestreams a set of n^ns processed maps — the 
annular transfer function maps F? ub (see Figure |4] for an 
example) — which show the response of the instrument 
at each pixel location to an annulus of constant temper- 
ature corresponding to bin b^. Finally, these maps are 
Hanning apodized, Fourier transformed, and represented 
as vectors F? ub of length n p i xe is- From these vectors, we 
construct an n p i xe i s x nbins matrix 

Bji h = (^ SUb )^ (5) 

the annular transfer function matrix, where Bjf h gives 
the response in the jth Fourier mode of the band sub- 



tracted map to an impulse in bin b{. 

Using this formalism, we can write a set of equations 
relating the Fourier transform of the Hanning apodized 
measured SPT map M| ub to the sky temperatures Tf. 

Mf b = B^ h Ti. (6) 

Since there are many more map pixels than radial bins, 
this is an over constrained system of linear equations, for 
which the least-squares solution is given by 

T = (B T WB)~ l B T WM, (7) 

where we have suppressed the superscript denoting band 
subtracted quantities and where W is the weight matrix. 
The weights are given by W = 1/(N + S CMB + S PS ), 
where TV is the detector and atmospheric noise and 
^cmb anc j gPS are ^ e CMB and dusty point source 
signals. These terms are all computed for the band sub- 
tracted map using the method described in Section l2.3.31 
Weighting the Fourier modes by their noise reduces 
the profile errors by 20%-30% compared with uniform 
weights. The covariance matrix for the profile is given by 
(B T WB)~ 1 , plus small additional terms to account for 
calibration and beam uncertainties. Figure [5] shows the 
importance of this procedure: the amplitudes and shapes 
of the profiles are badly misestimated if the beams and 
time-domain filtering are not taken into account. 

Since the SPT measures temperature differences, the 
cluster maps have arbitrary constant offsets. These man- 
ifest themselves as offsets in the profiles, and could bias 
the model fits if they were not taken into account. We 
include a constant offset as a parameter in our model fit- 
ting, and marginalize over this parameter in order to de- 
termine the quoted central decrements and angular scale 
radii. This marginalization ensures that our uncertainty 
in the zero point of the map is correctly accounted for in 
the model fits. The profiles shown in the Appendix are 
plotted with the offset fixed such that the temperature 
decrement in the outermost radial bin is zero. 

4.2. Model fitting 

We make use of the cluster profiles and their noise co- 
variance matrices to perform maximum likelihood fits to 
proposed models of the ICM. Our likelihood function is 
given by 

C(T\M) = 1 = exp , (8) 

V(2^)Met(AU ise ) V * J 

where n — 10 is the number of radial bins, T is the 
measured profile, M is the proposed model, and \ 2 is 
given by 

X 2 = (M - T) T (N noise )-\M - T). (9) 

We will consider two models: the /3-model and the 
GNFW model. 

The /?— model for thermal SZ cluster signals is moti- 
vated by the widely used /3— model parameterization of 
the three-dimensional electron number density: 

/ r 2 \ -3/3/2 

n e (r) =n eo I 1 + ^— ) , (10) 

V ^"core / 

where n e is the electron number density, n eo is the num- 
ber density at the cluster center, r is the radius from the 
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Fig. 4. — Bin 63 (left) and its 150 GHz annular transfer function map (right). Units are percent of maximum. The annular transfer 
function map represents the response of the SPT to a corresponding annulus of constant temperature, taking into account the beam and 
the time-domain filtering. 
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Fig. 5. — AS 1063 profile with and without beam and time- 
domain filtering correction. The uncorrected profile is computed 
by averaging the pixels in the band subtracted map within the 
radial bins, and the corrected profile is computed as described in 
Section O 

cluster center, r core is the core radius of the gas distribu- 
tion, and /3 is the power law index at large radii. Using 
this density function and assuming isothermality, Equa- 
tion [2] can be integrated along the line of sight through 
the cluster to give a simple analytic formula for the SZ 
decrement: 



AT SZ = AT 1 



2 



(l-3/3)/2 



(11) 



where 6 is the angular distance from the cluster cen- 
ter, given by = v/Da where Da is the angular diam- 
eter distance (computed assuming that (^m,^a^) = 
(0.3,0.7,0.7)), AXb is the central temperature decre- 
ment, and # CO re is the angle corresponding to r core . Al- 
though we do not necessarily expect the clusters in our 
sample to be isothermal, we nevertheless adopt Equa- 
tion[TT]as a convenient parameterization of the SZ signal. 

Two of the /3— model parameters, j3 and r core , are 
highly degenerate and diffi cult to constrain with SZ data 
alone (jGrego et al.ll200lD . We fix f3 at 0.86, the overall 
best-fit value for our sample (see Section [5^2]) . Typical X- 



ray clust er analyses yield val ues of f3 between about 0.6 
and 0.8 (|LaRoque et al.ll2006l and citations therein), but 
while the X-ray surfac e brightness is relat ively insensitive 
to ICM temperature (|Mohr et alJ Il999h , the SZ signal 
is proportional to the integrated pressure of the ICM. 
Therefore, any radial trend in temperature would lead 
to systematic differences between X-ray and SZ profiles. 
Our finding of f3 ~ 0.86 suggests that the ICM temper- 
ature tends to decrease with increasing radius, in qual- 
itative agreement with d irect measureme nts of cluster 
temperature profiles (e.g. JPratt et al.ll2QQ7h and w i th hy - 
drodynamical cluster simulations. Hal lman et al.l (|2007l ) 
find that j3 values for SZ profiles are systematically higher 
than X-ray /3 values by a factor of 1.21, consistent with 
the SPT results. 

We also fit to the GNFW model, which parameterizes 
the three-dimensional electron pressure as 



Pe(r) 



(r/r s )7n(l + (r/r s )-n)(^-T»)/«»' 



(12) 



where P eo is the pressure at the cluster center; r is the 
radius from the cluster center; r s = r 5 oo/c5oo is a scaling 
radius set by r 5 oo, the radius within which the average 
cluster density falls to 500 times the critical density of 
the universe at the cluster redshift, and C500, a parameter 
characterizing the gas concentration; and a n , /3 n , and j n 
set the slope at intermediate (r ~ r s ), large (r > r s ), and 
small (r < r s ) radii, respectively. We fix (a n , /? n > 7n> C500) 
to (1.0, 5.5, 0.5, 1.0), the parameters found to be the best 
fit to our stacked cluster profile (see Section [Q]) . These 
values are very similar to those found bv lN07l when fitting 
to Chandra data, (a n , /? n , 7 n , C500) = (0.9, 5.0, 0.4, 1.3) 27 . 
As the SZ effect is directly proportional to integrated 
pressure, we calculate ATsz given the pressure profile 
by numerically integrating the above function along the 
line of sight. As before, we define S = t s /Da- Note 
that fitting for r s with C500 fixed is equivalent to fitting 
for rsoo. However, because the quality of the GNFW 



27 These values are taken from Mroczkowski et al. (2009), and 
are updated from the results in 1N07I . 
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fit varies across our sample, we will later make use of a 
scaling relation to estimate rsoo- 

Both models thus contain two free parameters: a cen- 
tral temperature decrement, and a scaling radius (# C ore 
for the /?— model and 6 S for the GNFW model). As dis- 
cussed in Section l4TT| an offset parameter is also included. 

4.3. The effect of CMB and atmospheric noise 

The model fitting technique described in Section 14.21 
makes use of the full profile noise covariance matrices, 
including off-diagonal terms that describe correlations. 
This is important because the noise in the SPT maps can 
introduce spurious structure in the profiles that is highly 
correlated between bins. Most of these correlations arise 
from primary CMB anisotropics and atmospheric noise, 
the combination of which dominates the uncertainties in 
the profiles. 

The large angular scale confusion due to primary 
CMB anisotropics represents a fundamental limitation 
for single-band measurements of SZ radial profiles. The 
SPT 220 GHz data allow us to reduce the level of con- 
fusion; however, the 220 GHz maps have significantly 
more atmospheric noise than do the 150 GHz maps. 
Timestream filtering can help to mitigate atmospheric 
noise, but filtering too heavily will remove CMB signal. 
We chose our level of time-domian filtering by balancing 
these two effects, in order to maximize the CMB signal 
to noise at 220 GHz. 

The matched spatial filter described in Section 12.3.31 
is constructed to downweight noise-contaminated modes 
in the 220 GHz maps, and thus to increase the effective- 
ness of the CMB subtraction. Other spatial filters can 
be applied to the 220 GHz maps in order to change the 
balance between CMB and atmospheric confusion in the 
band subtracted map. We constructed several alterna- 
tive filters which improved the CMB subtraction at the 
expense of introducing atmospheric noise. In all cases, 
the resulting profiles and fits were found to be consis- 
tent but the errors were increased: by construction, the 
matched filter yields the best signal to noise in the band 
subtracted maps. Optimal band subtraction reduces the 
large angular scale confusion by a factor of 2 or more, 
depending on the depth of the 220 GHz map. 

Nevertheless, large-scale correlated noise due to at- 
mospheric fluctuations and unsubtracted CMB remains 
present in the band subtracted maps. Figure [6] shows 
the effect of this confusion on one band subtracted pro- 
file with particularly low 220 GHz signal to noise. While 
this cluster is an extreme example, the same issue per- 
sists across the entire sample. As in earlier SZ works, 
we find that model fits are in some respects more ro- 
bust indicators of global cluster properties than are the 
measured profiles, since the choice of a realistic model im- 
poses sensible constraints on less well constrained modes. 
We therefore make extensive use of model fitting when 
interpreting the SPT profiles and integrated Componiza- 
tions. 

5. RESULTS 

5.1. Profiles and model fits 

The projected radial profiles of the band subtracted 
maps can be found in the Appendix. Dashed vertical 
lines in the profile plots denote estimates of r$oo deter- 
mined from the temperature-radius scaling relation in 
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Fig. 6. — Measured profile (black) and best-fit /3-model (red) 
for RXCJ2031. 8-4037. The noise level in the 220 GHz map for this 
cluster is among the highest in the sample, so the noise correlations 
between bins are large. Shown in gray are the results of varying 
the first component of the measured profile by 1, 2, and 3a in the 
eigenbasis of the noise covariance matrix. 

IVikhlinin etHI (|2006[ ). The error bars on the profiles 
are determined using Monte Carlo methods, and as dis- 
cussed in Section 14. 3| are strongly correlated between 
bins. 

The best-fit central decrement and scaling radius for 
each cluster are listed in Table [3j along with the mini- 
mum x 2 an d the logarithmic ratio of the marginal likeli- 
hoods of the GNFW and f3— models: 

i? = log 10 (^J C GNFW d6 s dAT / J £pd0 core dAT^J 

(13) 

where £gnfw is the likelihood function for the GNFW 
model and Cp is the likelihood function for the (3— model. 
This quantity can be interpreted as a likelihood ratio, 
with a value of 10 corresponding to a 10:1 preference for 
the GNFW model. For the recent mergers 1ES 0657-56 
and A 2744, neither model is a good fit, though the /?- 
model is strongly favored. The f3— model yields a larger 
marginal likelihood for over half of the sample, but the 
SPT data do not definitively favor one model over the 
other. 

5.2. Stacked profiles 

Each individual cluster profile contains too little in- 
formation to place a tight constraint on the pressure be- 
havior at large radius. We proceed by scaling the profiles 
to rsoo , as det ermined from the temp erature-radius scal- 
ing relation in IVikhlinin et al.l (|2006f ) , and analyzing the 
data set as a whole — by constructing a stacked profile, 
and by fitting models to the sample. The stacked profile, 
shown in Figure [7J is determined by averaging the scaled 
profiles weighted by their noise. Each radial bin in the 
stacked profile is assigned a radius and a temperature 
from the average of the scaled profiles, and an uncer- 
tainty is assigned by the standard error on the mean. 

We fit models to the sample by evaluating the likeli- 
hood function for each individual cluster at each point in 
a grid of parameter space, and then by multiplying the 
cluster likelihoods together to produce a joint likelihood. 
This allows us to determine average values for the struc- 
ture parameters in our two models. For the /3- model, 
we allow f3 and r core /r^oo to vary. We allow the cen- 
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TABLE 3 
Best-fit Model Parameters 





Isothermal (3 — model 


GNFW Model 




ID 


AT ( M K CMR ) 


#core (') 


AT ( M K CMR ) 


S (') 


R 




PTE 


Primary sample 


A 2744 


-549 ± 27 


1.51 ± 0.06 


-670 ± 33 


7.87 ± 0.51 


-50.79 


19.50 


0.01 


RXCJ0217. 2-5244 


-214 ± 29 


0.64 ± 0.13 


-262 ± 34 


3.12 ± 0.63 


-0.84 


5.24 


0.63 


RXCJ0232. 2-4420 


-510 ± 26 


0.99 ± 0.05 


-626 ± 32 


5.01 ± 0.31 


17.72 


9.68 


0.21 


AS 0520 


-217 ± 13 


1.22 ± 0.12 


-268 ± 16 


6.80 ± 0.86 


16.33 


9.16 


0.24 


RXCJ0528. 9-3927 


-335 ± 20 


1.78 ± 0.14 


-412 ± 25 


10.03 ± 1.24 


5.53 


7.46 


0.38 


AS 0592 


-529 ± 31 


0.79 ± 0.06 


-658 ± 40 


3.77 ± 0.29 


-3.51 


3.00 


0.89 


A 3404 


-472 ± 29 


0.95 ± 0.08 


-588 ± 38 


4.52 ± 0.43 


1.20 


4.09 


0.77 


1ES 0657-56 


-932 ± 43 


1.47 ± 0.03 


-1133 ± 52 


7.56 ± 0.22 


-270.59 


65.61 


0.00 


RXCJ2031. 8-4037 


-511 ± 42 


0.65 ± 0.07 


-616 ± 48 


3.24 ± 0.35 


-3.76 


11.66 


0.11 


A 3888 


-516 ± 26 


1.35 ± 0.07 


-634 ± 32 


6.79 ± 0.49 


-28.25 


8.64 


0.28 


AS 1063 


-1062 ± 50 


0.86 ± 0.02 


-1321 ± 62 


4.14 ± 0.12 


-2.54 


28.61 


0.00 


Supplemental sample 


RXCJ0336. 3-4037 


-306 ± 30 


0.70 ± 0.10 


-388 ± 40 


3.20 ± 0.46 


-5.11 


28.33 


0.00 


RXCJ0532. 9-3701 


-449 ± 24 


1.00 ± 0.07 


-543 ± 29 


5.66 ± 0.52 


21.25 


35.43 


0.00 


M ACS J0553. 4-3342 


-680 ± 34 


0.96 ± 0.05 


-849 ± 44 


4.54 ± 0.26 


-18.63 


28.33 


0.00 


A 3856 


-256 ± 26 


0.71 ± 0.11 


-313 ± 31 


3.49 ± 0.55 


-0.01 


2.17 


0.95 



Note. — The /3-model fit is performed with (3 fixed to 0.86, and the GNFW fit is performed with the parameters (ain,/3n,7n, C500) fixed 
to (1.0, 5.5, 0.5, 1.0). The sixth column quantifies the relative goodness of fit in terms of the marginal likelihood ratio (see the text). Positive 
values indicate that the GNFW model is preferred, and negative values indicate that the (3— model is preferred. The seventh column is the 
X 2 for the preferred model, and the eighth is the probability to exceed this value of % 2 for seven degrees of freedom. 



tral decrement and overall constant offset of each profile 
to vary as well, and we marginalize over these param- 
eters. The cluster center position is held fixed, leaving 
six remaining degrees of freedom per cluster. We find 
that for the full sample, the marginalized values of j3 
and r cor e/^500 are 0.86 ±0.09 and 0.20 ±0.05. Excluding 
the supplemental sample from the fit does not appre- 
ciably change the results. The minimum x 2 for the fit 
is 338 for 90 degrees of freedom, indicating a relatively 
poor fit. If f3 is fixed to a typical X-ray value of 2/3, 
the minimum x 2 is 362. If 1ES 0657-56 is excluded from 
the fit, the minimum reduced x 2 changes from 338/90 = 
3.76 to 296/84 = 3.52, and the marginalized values of 
(/?, r core /r 50 o) change to (0.78 ± 0.08, 0.17 ± 0.04). 

For the generalized NFW model, we fit for the struc- 
ture parameters a n , /3 n , 7 n , and C500, again fixing the 
center positions and marginalizing over the central decre- 
ments of all of the clusters. The likelihood- weighted 
mean values are (a n , /3 n , 7 n , C500) = (1-0, 5.5, 0.5, 1.0) and 
the marginalized values are (1.2 ± 0.5,5.0 ± 2.0,0.7 ± 
0.3,1.3 ± 0.5). The mini mum y 2 for the fit is 468, 
compared to 472 for the lN07l parameter values of 
(a n ,/3n,7n,c 5 oo) = (0.9,5.0,0.4, 1.3)— the change in x 2 
is small due to strong degeneracies between the param- 
eters. If 1ES 0657-56 is excluded from the GNFW 
fit, the marginalized parameter values do not change 
appreciably, but the minimum reduced x 2 drops from 
468/75 = 6.24 to 318/70 = 4.54. 

In Figure El the SPT maximum likelihood GNFW 
model with (a n , /3 n , 7 n , C5oo) = (l-0, 5.5, 0.5, 1.0) is com- 
pared to the (0.9,5.0,0.4,1.3) model foun d bv lN07l and 
the (1 .0 5,5.49 ,0.31,1.18) model found bv lArnaud et al.l 
(200 91). IN07I use a samp le of nearby relaxed clusters, 
while lArnaud et al.l ([2009) use the REXCESS represen- 
tative sample of z < 0.2 clusters. Two systems — A 3888 
and A 3856— are in both the SPT and REXCESS sam- 
ples. On average, however, the samples are quite dif- 
ferent, with the SPT clusters tending to be hotter and 
more distant. The similarity between the resulting pro- 
files is encouraging for future applications of self-similar 
pressure models. 
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Fig. 7. — Stacked radial profile and maximum likelihood (3- 
model and GNFW model fits. The plot extends to 2r5oo, which 
is approximately equal to the average virial radius for the sample. 
The best-fit model profiles are computed from simulated processed 
maps i n the same way that the cluster profiles are computed (see 
Section |4.1|) . and are shown as a continuous line only for the pur- 
pose of visualization. 

The SPT best-fit /3-model and GNFW model are 
shown in Figure [7J Since the models are very similar, 
the fits are of roughly equal quality. The logarithmic 
ratio of marginal likelihoods (see Equation [T3|) for the 
GNFW model compared to the /3-model is R = —322, 
indicating a preference for the latter. This preference 
is due primarily to 1ES 0657-56, and to a lesser extent 
to A 2744. Both of these clusters are major mergers. 
Excluding the former brings the marginal likelihood ra- 
tio to R = —51, and excluding both brings the ratio to 
R = —0.6. Stronger model discrimination will require 
even better multi-band measurements to large angular 
scales, where atmospheric noise and astrophysical con- 
tamination limit the constraining power of the data. 

5.3. Integrated y parameter and scaling relations 

The integrated Compton-?/ parameter >sz is expected 
to serve as a good proxy for mass, and to be related 
to other physical properties of clusters by scaling rela- 
tions. Ysz can be computed either using the model fits 
discussed in Section [4721 or using the model-independent 
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Fig. 8. — GNFW mode ls with be st-fit slope and concentration 
parameters from N07, Arnaud et al. (2009), and this work. Note 
that the samples in each work differ in characteristic redshift, mor- 
phology, and temperature. These results are indicative of the con- 
sistency between X-ray and SZ pressure measurements. 

radial profiles. Since the large angular scale noise con- 
tributes strongly to the model-independent estimate, we 
consider the model-dependent results to be slightly more 
robust. The average ratio across the sample of Ysz com- 
puted from the f3- model fit, to Ysz computed from the 
GNFW model fit, is 1.06 ± 0.08 at r 25 oo and 1.19 ± 0.20 
at rsoo- This is due to the fact that a /3 = 0.86 model 
falls off more gradually at intermediate radius than does 
the GNFW model with our best-fit structure parameters. 
Since the /3-model is slightly preferred over the GNFW 
model across our sample, we use the f3- model fits to com- 
pute our model-dependent results; the relations derived 
from the GNFW values are consistent to within la. 

Both the model-dependent and the model-independent 
estimates of the integrated Compton-?/ parameter are 
given in Table HJ Uncertainties are determined via Monte 
Carlo methods, and the apertures are again determined 
from the temperature-radius scaling relation. Since the 
slope of the temperature-radius scaling relation is rela- 
tively shallow, the aperture radii vary only by ~ 10% 
across our sample. Although our integration radii are 
slightly different than those used in the X-ray analyses, 
we find that this choice makes a negligible difference in 
the integrated Y parameters — even integrating out to a 
fixed radius of 1 Mpc does not significantly increase the 
scatter in the scaling relations. 

The /3-model and model-independent estimates at r25oo 
are in relatively good agreement, with a ratio of model- 
dependent to model-independent Ysz of 0.99 ± 0.22. 
Significant scatter is seen at rsoo, where the ratio is 
1.30 ±0.56. The model-independent SZ signal is less well 
constrained at larger radius due to its lower amplitude 
and to the large angular scale confusion discussed in Sec- 
tion [4731 The model-dependent estimates are less affected 
by confusion, as noisier modes are downweighted in the 
model fitting process. We therefore use the /3-model es- 
timates to explore scaling relations with X-ray values, 
w hich are also der ive d in a model-dependent manner. 

lHalverson et al.l ([2009) have recently observed 
1ES 0657-56 with the APEX-SZ experiment, and 
have fit an ellipsoidal /3-model to their 150 GHz SZ 
map. Since we assume azimuthal symmetry and fit 
models to the radial profiles rather than the maps, 
we cannot perform a directly comparable fit. How- 



ever, we can compare our respective estimates of 
Ysz. The APEX-SZ best-fit ellipsoidal model yields 
(*sz,2500,*sz,50o) = (30 ± 5, 66 ± 18) x lO" 11 sr, 
which is consistent at the ~ la level with both the 
model-dependent and model-independent SPT results. 

The /3-model estimates of Ysz, 500 are plotted ver- 
sus X-ray estimates of gas mass and Yx in Figure [9l 
Yx = M ga , a fcT R is an X-ray observable analogous to Ysz 
([Kravtsov et al.l l2006). Only clusters with published val- 
ues of gas mass and electron temperature from the refer- 
ences in Table [T] are included in the plots. Several clus- 
ters have multiple published X-ray estimates of electron 
temperature and gas mass. For our plot s in Figure [H 
we chose to adopt the I Zhang et al.l ([2006, 2008) results 
where available, since these samples had the largest over- 
lap with our own. For different choices of published X- 
ray values, we find scaling relation parameters that are 
consistent at the 2a level. 

In order to compare Yx and Ysz, we convert the former 
from a spherically defined quantity into a cylindrical pro- 
j ection analogous to Ysz as described in lBonamente et ail 
(2008). This method requires an estimate of the ratio of 
the cluster pressure within a cylindrical and spherical vol- 
ume, which we calculate assuming the best-fit SZ profile 
for each cluster. Since this ratio has large uncertainties in 
the model-independent case, we report scaling relations 
only for our model-dependent Ysz estimates. Temper- 
ature structure in cl usters typically causes Yx to be bi- 
ased high by - 10% (|Arnaud et al Jl2QQ9t IVikhlinin et al.l 
2006); however, we consider this bias to be subdominant 
to the other uncertainties in this comparison. 

Ysz is also expected to scale with electron temperature 
and total mass. However, the reported errors in the X- 
ray estimates of electron temperature vary widely across 
our sample, and different analyses of the same clusters 
often produce inconsistent results. The X-ray total mass 
estimates are determined under the assumption of hy- 
drostatic equilibrium, and the errors are quite large. For 
this reason, we do not include the total mass and electron 
temperature scaling relations in this work. In a future 
paper, we will pursue a joint SZ and X-ray analysis which 
can more appropriately incorporate these uncertainties. 

The Ysz -gas mass plot is in a form readily comparabl e 
to the scaling relation given in lBonamente et all ((2008): 

Y SZ D\ cx /-J/X/s 3 ^) 273 , (14) 
where Da is the angular diameter distance, / gas is the 
gas mass fraction (which we assume to be constant), E(z) 
is given by ^/QmO- + z) 3 + ^a, and M gas is the cluster 
gas mass. We fit for the gas mass scaling relation in log- 
l og sp ace using the method described in iMarrone et al.l 
(|2009[ ) , assuming a scaling in the form Y = a + j3X and 
allowing for intrinsic scatter. We find a = —5.73 ± 0.43 
and /3 = 2.12 ± 0.45. The best-fit intrinsic scatter in this 
scaling relation, in terms of a percentage change in M gas 
at fixed Ysz-, is 14 ± 10%. If we substitute Ysz, 2500 f° r 
^5Z,500j we find a similar scaling relation with M gaS)5 oo' 
a = -5.92 ± 0.41, f3 = 1.97 ± 0.44, and a scatter of 
15% ±11%. 

The gas masses of two clusters in our sample, 
1ES 0657-56 and RXC J0232. 2-4420, ar e estim ated both 
from Chandra data in Maughan e t al.l (|2008l ) and from 
XMM data in I Zhang et al.l (|2006l ). In both cases, the 
iMaughan etHI ([2008) estimates (2.31±o.Sl x 10 14 M o for 
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TABLE 4 
Integrated Compton-2/ parameters. 



ID 


^2500 


^500 


f3— model 


GNFW Model 


Model-independent 








^SZ,2500 


^SZ,500 


^SZ,2500 


^SZ,500 


^SZ,2500 


^SZ,500 








Primary samp] 


e 








A 97/1/1 


^ 
U.DO 


1 /IO 


20.3 ±0.7 


49.2 ± 2.3 


19.8 ± 1.0 


47.8 ± 3.8 


19.0 ± 1.4 


38.1 ±6.8 




fi1 


l.oO 


3.3=1=0.4 


6.3 ±1.0 


3.1 ±0.4 


5.2 ±1.0 


3.6 ±0.9 


4.8 ±2.7 




^9 


1 17 


11.2 ±0.4 


24.8 ± 1.1 


10.8 ±0.4 


23.1 ± 1.5 


12.9 ±0.9 


26.7 ±4.3 


AS 0520 


0.54 


1.21 


5.7 ±0.4 


13.7 ± 1.4 


5.9 ±0.6 


14.4 ±2.1 


4.7 ±0.8 


7.2 ±2.7 


RXCJ0528.9-3927 


0.53 


1.19 


11.1 ±0.7 


31.0 ±2.9 


11.5 ±0.9 


33.0 ±4.3 


9.8 ±0.9 


22.1 ±4.4 


AS 0592 


0.58 


1.30 


12.9 ±0.6 


24.7 ± 1.6 


11.6 ±0.7 


18.7 ± 1.7 


14.7 ±2.0 


26.0 ±9.8 


A 3404 


0.60 


1.35 


18.2 ± 1.3 


34.6 ±3.1 


15.9 ± 1.5 


24.5 ±3.1 


13.5 ±3.0 


14.0 ±5.4 


1ES 0657-56 


0.65 


1.44 


36.2 ±0.6 


84.9 ± 1.8 


35.1 ±0.7 


81.0 ±2.4 


30.0 ± 1.0 


45.9 ±4.7 


RXCJ2031. 8-4037 


0.64 


1.43 


8.4 ±0.5 


16.1 ± 1.4 


7.7 ±0.5 


12.8 ± 1.3 


14.3 ± 1.6 


37.0 ±8.0 


A 3888 


0.59 


1.33 


30.1 ± 1.5 


60.8 ±3.7 


28.7 ±2.0 


51.9 ±5.2 


24.6 ±3.8 


36.1 ± 14.7 


AS 1063 


0.65 


1.44 


23.3 ±0.4 


48.3 ± 1.1 


21.5 ±0.5 


39.7 ± 1.5 


25.5 ± 1.1 


44.5 ± 5.3 


Supplemental sample 


RXCJ0336.3-4037 


0.56 


1.26 


7.5 ±0.7 


13.7 ± 1.7 


6.2 ±0.8 


8.9 ± 1.7 


11.0 ±2.5 


27.8 ± 11.8 


RXCJ0532.9-3701 


0.62 


1.38 


12.8 ±0.7 


26.6 ± 1.9 


13.6 ± 1.0 


27.6 ±3.2 


18.1 ± 1.7 


42.4 ±8.0 


MACSJ0553.4-3342 


0.69 


1.52 


15.8 ±0.5 


34.3 ± 1.5 


14.5 ±0.6 


28.8 ±2.0 


14.7 ± 1.4 


20.5 ±5.1 


A 3856 


0.55 


1.23 


7.3 ±0.8 


13.1 ± 1.7 


6.2 ±0.9 


8.4 ± 1.7 


7.7 ±2.9 


9.8 ±6.4 



Note. — Units for Ysz are 10 11 sr. Both model-dependent and model-independent estimates are provided (see 
the text). r25oo and rsoo are determined from the temperature-radius scaling relation. 



1ES 0657-56 and 1.201 ^ x 1O 14 M for RX CJ0232.2- 
4420) are larger than the Zha ng et al.l (j2QQ6h results we 
have adopted. If we instead adopt the larger values, we 
find /? = 1.85 ± 0.43, consistent with the self-similar ex- 
pectation. 

The Ysz-Yx relation is expected to have a small scat- 
ter and to have a slope of 1.0. Our best-fit power law 
parameters are a = 2.42 ± 0.87 and j3 = 1.56 ± 0.21, 
with an intrinsic scatter of 12% ±12% in Yx- Substi- 
tuting Fsz,2500 for Y S z,500 gives a scaling relation of 
a = 1.69±0.78, /? = 1.46±0.20, and a scatter of 9% ±9%. 
While these slopes are inconsistent with the expected 
value of 1.0, the discrepancy is largely due to the three 
lowest signal to noise data points: A 3856, AS 0520, and 
A 3404. In the gas mass scaling relation, we found that 
adopting independent published values of gas mass intro- 
duced significant variations in the slope and offset. Since 
Yx is computed from the same M gas estimates, compara- 
ble systematic uncertainties and offsets are expected. We 
therefore hesitate to draw conclusions before performing 
a more uniform X-ray analysis. 

We consider these F^-gas mass and Ysz~Yx scaling 
relations to be reasonably consistent with the self-similar 
values, and we find no measurable scatter within the 
limits of our statistical and s ystem atic errors. Previous 
studies (e.g.. iKraytsov et al.ll20q6t lArnaud et al.l 120071 : 
iNagai et al.ll200rt IVikhlinin et al.ll2QQ9h have found that 
M gas and Yx scale with the total cluster mass with low 
scatter, and these results indicate that Ysz as measured 
with the SPT should behave similarly. Future joint anal- 
yses of X-ray and SZ data will further constrain these 
scaling relations, and will provide better limits on any 
departure from self-similarity. 

6. CONCLUSIONS 

In this paper, we have summarized SZ observations 
of 15 bright galaxy clusters taken with the SPT. The 



SZ signals were measured at 150 GHz, and concurrent 
220 GHz measurements were used to remove astrophys- 
ical contamination. Radial profiles were computed for 
each cluster using a technique that accounts for the ef- 
fects of the beam and the time-domain filtering, and si- 
multaneously characterizes the correlated errors due to 
detector and atmospheric noise. 

The profiles were fit both to a j3— model and to a gen- 
eralized NFW model (|N07l ). By scaling and stacking 
the SPT profiles, and allowing the structure parameters 
of the models to vary, we obtained self-similar pressure 
profiles that agree closely with previous X-ray and SZ 
results. For the j3— model, we found j3 — 0.86 ± 0.09 
and r C ore/r500 = 0.20 ± 0.05. For the GNFW model, 
we found (a n , /3 n , 7m C5qq ) = (1 .0, 5. 5, 0.5, 1.0), similar t o 
the parameters found bv lN07l and I Arnaud et aD ((2009). 
The SPT data do not strongly prefer one model over the 
other. Our measurements show no significant difference 
between the stacked SZ profile and the models out to 
beyond rsoo- This is the first SZ-measured constraint on 
cluster pressure profiles at such large radii. 

We also estimated the integrated Compton-?/ parame- 
ter Ysz f° r each cluster using both model-dependent and 
model-independent techniques. The scaling relations be- 
tween Ysz and the X-ray-determined M gas and Yx were 
found to be reasonably consistent with the self-similar 
values. We find an intrinsic scatter in the Yx~Ysz scaling 
relation of 9% ±9% at r 25 oo and 12% ±12% at r 50 o. Pre- 
vious studies have indicated that gas mass and Yx scale 
with low scatter to the total cluster mass. These results 
indicate that Ysz as measured with the SPT should be- 
have similarly, which is encouraging for the ongoing SPT 
cluster SZ survey. 

The data presented in this work demonstrate the util- 
ity of SZ measurements for characterizing the ICM out to 
large radii. In future works, we will expand the SPT sam- 
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Fig. 9. — Comparisons between SZ and X-ray measurements within r^oo: Ys^-Mg as (left) and Ysz~Yx (right). X-ray values are taken 
from the references in Tabled and clusters without published values in these sources are omitted. SZ values are taken from the /3-model 
estimates in Table [4[ In the right panel, Cx — C{<tt /m e c 2 )/ fi e Tn p where C is a constant accounting for spherical versus cylindrical 
integration (Bonamente et al. 2008). The dashed lines represent the best-fit power laws, and the shaded area represents the la allowed 
region. The dotted lines represent the best bit powers law with the slopes fixed to the expected self-similar values of 5/3 (left) and 1 (right). 



pie of X-ray luminous clusters, and will include 95 GHz 
data in addition to 150 and 220 GHz. The addition of a 
third frequency band will allow us to better remove astro- 
physical backgrounds such as the CMB, thus improving 
our estimates of thermal SZ signals at large radii. Future 
analyses will also combine SZ measurements with X-ray 
cluster observations in order to separately estimate the 
temperature and density of the ICM, and to determine 
the cluster gas mass fractions. 

The South Pole Telescope is supported by the National 
Science Foundation through grants ANT-0638937 and 
ANT-0130612. Partial support is also provided by the 
NSF Physics Frontier Center grant PHY-01 14422 to the 



Kavli Institute of Cosmological Physics at the University 
of Chicago, the Kavli Foundation and the Gordon and 
Betty Moore Foundation. The McGill group acknowl- 
edges funding from the National Sciences and Engineer- 
ing Research Council of Canada, the Quebec Fonds de 
recherche sur la nature et les technologies, and the Cana- 
dian Institute for Advanced Research. The following 
individuals acknowledge additional support: B. Benson 
and K. Schaffer from a KICP Fellowship; J. McMahon 
from a Fermi Fellowship; Z. Staniszewski from a GAAN 
Fellowship; A. T. Lee from the Miller Institute for Basic 
Research in Science, University of California, Berkeley; 
and N. W. Halverson from an Alfred P. Sloan Research 
Fellowship. 



REFERENCES 



Arnaud, M., Pointecouteau, E., & Pratt, G. W. 2007, A&A, 474, 
L37 

Arnaud, M., Pratt, G. W., Piffaretti, R., Boehringer, H., Croston, 

J. H., & Pointecouteau, E. 2009, ArXiv e-prints 
Benson, B. A., Church, S. E., Ade, P. A. R., Bock, J. J., Ganga, 

K. M., Henson, C. N., & Thompson, K. L. 2004, ApJ, 617, 829 
Birkinshaw, M. 1999, Physics Reports, 310, 97 
Bohringer, H., Schuecker, P., Guzzo, L., Collins, C. A., Voges, W., 

Cruddace, R. G., Ortiz-Gil, A., Chincarini, G., De Grandi, S., 

Edge, A. C, MacGillivray, H. T., Neumann, D. M., Schindler, 

S., & Shaver, P. 2004, A&A, 425, 367 
Bonamente, M., Joy, M., LaRoque, S. J., Carlstrom, J. E., Nagai, 

D., & Marrone, D. P. 2008, ApJ, 675, 106 
Borys, C, Chapman, S., Halpern, M., & Scott, D. 2003, MNRAS, 

344, 385 

Bussmann, R. S., Holzapfel, W. L., & Kuo, C. L. 2005, ApJ, 622, 
1343 

Calabretta, M. R. & Greisen, E. W. 2002, A&A, 395, 1077 
Carlstrom, J. E., Ade, P. A. R., Aird, K. A., Benson, B. A., Bleem, 
L. E., Busetti, S., Chang, C. L., Chauvin, E., Cho, H., Crawford, 
T. M., Crites, A. T., Dobbs, M. A., Halverson, N. W., Heimsath, 
S., Holzapfel, W. L., Hrubes, J. D., Joy, M., Keisler, R., Lanting, 
T. M., Lee, A. T., Leitch, E. M., Leong, J., Lu, W., Lueker, M., 
McMahon, J. J., Mehl, J., Meyer, S. S., Mohr, J. J., Montroy, 
T. E., Padin, S., Plagge, T., Pryke, C, Ruhl, J. E., Schaffer, 
K. K., Schwan, D., Shirokoff, E., Spieler, H. G., Staniszewski, 
Z., Stark, A. A., & Vieira, K. V. J. D. 2009, submitted to PASP, 
arXiv:0907.4445 

Carlstrom, J. E., Holder, G. P., & Reese, E. D. 2002, ARA&A, 40, 
643 



Cavagnolo, K., Donahue, M., Voit, G., & Sun, M. 2008, The 

Astrophysical Journal, 682, 821 
Cavaliere, A. & Fusco-Femiano, R. 1976, A&A, 49, 137 
Cavaliere, A. & Fusco-Femiano, R. 1978, A&A, 70, 677 
Clowe, D., Bradac, M., Gonzalez, A. H., Markevitch, M., Randall, 

S. W., Jones, C, & Zaritsky, D. 2006, ApJ, 648, L109 
Colafrancesco, S., de Bernardis, P., Masi, S., Polenta, G., & Ullio, 

P. 2007, A&A, 467, LI 
da Silva, A. J. C. 2004, Ap&SS, 290, 167 

Finoguenov, A., Bohringer, H., & Zhang, Y.-Y. 2005, A&A, 442, 
827 

Grego, L., Carlstrom, J. E., Joy, M. K., Reese, E. D., Holder, G. P., 
Patel, S., Cooray, A. R., & Holzapfel, W. L. 2000, ApJ, 539, 39 

Grego, L., Carlstrom, J. E., Reese, E. D., Holder, G. P., Holzapfel, 
W. L., Joy, M. K., Mohr, J. J., & Patel, S. 2001, ApJ, 552, 2 

Haehnelt, M. G. & Tegmark, M. 1996, MNRAS, 279, 545+ 

Hall, N. R., Knox, L., Reichardt, C. L., Ade, P. A. R., Aird, K. A., 
Benson, B. A., Bleem, L. E., Carlstrom, J. E., Chang, C. L., 
Cho, H., Crawford, T. M., Crites, A. T., de Haan, T., Dobbs, 
M. A., George, E. M., Halverson, N. W., Holder, G. P., Holzapfel, 
W. L., Hrubes, J. D., Joy, M., Keisler, R., Lee, A. T., Leitch, 
E. M., Lueker, M., McMahon, J. J., Mehl, J., Meyer, S. S., Mohr, 
J. J., Montroy, T. E., Padin, S., Plagge, T., Pryke, C, Ruhl, 
J. E., Schaffer, K. K., Shaw, L., Shirokoff, E., Spieler, H. G., 
Staniszewski, Z., Stark, A. A., Switzer, E. R., Vanderlinde, K., 
Vieira, J. D., Williamson, R., & Zahn, O. 2009, submitted to 
ApJ, arXiv:0912.4315 

Hallman, E. J., Burns, J. O., Motl, P. M., & Norman, M. L. 2007, 
ApJ, 665, 911 



15 



Halverson, N. W., Lanting, T., Ade, P. A. R., Basu, K., Bender, 
A. N., Benson, B. A., Bertoldi, F., Cho, H., Chon, G., Clarke, J., 
Dobbs, M., Ferrusca, D., Glisten, R., Holzapfel, W. L., Kovacs, 

A. , Kennedy, J., Kermish, Z., Kneissl, R., Lee, A. T., Lueker, 
M., Mehl, J., Menten, K. M., Muders, D., Nord, M., Pacaud, F., 
Plagge, T., Reichardt, C., Richards, P. L., Schaaf, R., Schilke, 
P., Schuller, F., Schwan, D., Spieler, H., Tucker, C., Weiss, A., 
& Zahn, O. 2009, ApJ, 701, 42 

Herranz, D., Sanz, J. L., Barreiro, R. B., & Martinez- Gonzalez, E. 

2002a, ApJ, 580, 610 
Herranz, D., Sanz, J. L., Hobson, M. P., Barreiro, R. B., Diego, 

J. M., Martinez-Gonzalez, E., & Lasenby, A. N. 2002b, MNRAS, 

336, 1057 

Hincks, A. D., Acquaviva, V., Ade, P., Aguirre, P., Amiri, M., 
Appel, J. W., Barrientos, L. F., Battistelli, E. S., Bond, J. R., 
Brown, B., Burger, B., Chervenak, J., Das, S., Devlin, M. J., 
Dicker, S., Doriese, W. B., Dunkley, J., Diinner, R., Essinger- 
Hileman, T., Fisher, R. P., Fowler, J. W., Hajian, A., Halpern, 
M., Hasselneld, M., Hernandez-Monteagudo, C., Hilton, G. C., 
Hilton, M., Hlozek, R., Huffenberger, K., Hughes, D., Hughes, 
J. P., Infante, L., Irwin, K. D., Jimenez, R., Juin, J. B., Kaul, M., 
Klein, J., Kosowsky, A., Lau, J. M., Limon, M., Lin, Y., Lupton, 
R. H., Marriage, T., Marsden, D., Martocci, K., Mauskopf, P., 
Menanteau, F., Moodley, K., Moseley, H., Netterneld, C. B., 
Niemack, M. D., Nolta, M. R., Page, L. A., Parker, L., Partridge, 

B. , Quintana, H., Reid, B., Sehgal, N., Sievers, J., Spergel, D. N., 
Staggs, S. T., Stryzak, O., Swetz, D., Switzer, E., Thornton, R., 
Trac, H., Tucker, C., Verde, L., Warne, R., Wilson, G., Wollack, 
E., & Zhao, Y. 2009, ArXiv e-prints 

Hinshaw, G., Weiland, J. L., Hill, R. S., Odegard, N., Larson, D., 
Bennett, C. L., Dunkley, J., Gold, B., Greason, M. R., Jarosik, 
N., Komatsu, E., Nolta, M. R., Page, L., Spergel, D. N., Wollack, 

E. , Halpern, M., Kogut, A., Limon, M., Meyer, S. S., Tucker, 
G. S., & Wright, E. L. 2009, ApJS, 180, 225 

Hughes, J. P., Menanteau, F., Sehgal, N., Infante, L., &; Barrientos, 

F. 2009, in Bulletin of the American Astronomical Society, 
Vol. 41, Bulletin of the American Astronomical Society, 336 — h 

Itoh, N., Kawana, Y., Nozawa, S., & Kohyama, Y. 2000, ApJ 
Kravtsov, A. V., Vikhlinin, A., & Nagai, D. 2006, ApJ, 650, 128 
LaRoque, S. J., Bonamente, M., Carlstrom, J. E., Joy, M. K., 

Nagai, D., Reese, E. D., &; Dawson, K. S. 2006, ApJ, 652, 917 
Lima, M., Jain, B., &; Devlin, M. 2009, ArXiv e-prints 
Markevitch, M., Forman, W. R., Sarazin, C. L., & Vikhlinin, A. 

1998, ApJ, 503, 77 
Marrone, D. P., Smith, G. P., Richard, J., Joy, M., Bonamente, 
M., Hasler, N., Hamilton-Morris, V., Kneib, J., Culverhouse, T., 
Carlstrom, J. E., Greer, C, Hawkins, D., Hennessy, R., Lamb, 
J. W., Leitch, E. M., Loh, M., Miller, A., Mroczkowski, T., 
Muchovej, S., Pryke, C, Sharp, M. K., & Woody, D. 2009, ApJ, 
701, L114 

Maughan, B. J., Jones, C, Forman, W., &; Van Speybroeck, L. 

2008, ApJS, 174, 117 
Melin, J.-B., Bartlett, J. G., & Delabrouille, J. 2006, A&A, 459, 

341 

Mohr, J. J., Mathiesen, B., & Evrard, A. E. 1999, ApJ, 517, 627 
Motl, P. M., Hallman, E. J., Burns, J. O., & Norman, M. L. 2005, 
ApJ, 623, L63 

Mroczkowski, T., Bonamente, M., Carlstrom, J. E., Culverhouse, 
T. L., Greer, C, Hawkins, D., Hennessy, R., Joy, M., Lamb, 
J. W., Leitch, E. M., Loh, M., Maughan, B., Marrone, D. P., 
Miller, A., Muchovej, S., Nagai, D., Pryke, C, Sharp, M., & 
Woody, D. 2009, ApJ, 694, 1034 

Nagai, D. 2006, ApJ, 650, 538 

Nagai, D., Kravtsov, A. V., &; Vikhlinin, A. 2007, ApJ, 668, 1 
Nolta, M. R., Dunkley, J., Hill, R. S., Hinshaw, G., Komatsu, E., 
Larson, D., Page, L., Spergel, D. N., Bennett, C. L., Gold, B., 
Jarosik, N., Odegard, N., Weiland, J. L., Wollack, E., Halpern, 
M., Kogut, A., Limon, M., Meyer, S. S., Tucker, G. S., & Wright, 
E. L. 2009, ApJS, 180, 296 
Padin, S., Staniszewski, Z., Keisler, R., Joy, M., Stark, A. A., Ade, 
P. A. R., Aird, K. A., Benson, B. A., Bleem, L. E., Carlstrom, 
J. E., Chang, C. L., Crawford, T. M., Crites, A. T., Dobbs, M. A., 
Halverson, N. W., Heimsath, S., Hills, R. E., Holzapfel, W. L., 
Lawrie, C, Lee, A. T., Leitch, E. M., Leong, J., Lu, W., Lueker, 
M., McMahon, J. J., Meyer, S. S., Mohr, J. J., Montroy, T. E., 
Plagge, T., Pryke, C, Ruhl, J. E., Schaffer, K. K., Shirokoff, E., 
Spieler, H. G., & Vieira, J. D. 2008, Appl. Opt., 47, 4418 



Piffaretti, R., Jetzer, P., Kaastra, J. S., & Tamura, T. 2005, A&A, 
433, 101 

Pratt, G. W., Bohringer, H., Croston, J. H., Arnaud, M., Borgani, 
S., Finoguenov, A., &; Temple, R. F. 2007, A&A, 461, 71 

Reese, E. D., Carlstrom, J. E., Joy, M., Mohr, J. J., Grego, L., & 
Holzapfel, W. L. 2002, ApJ, 581, 53 

Reichardt, C. L., Ade, P. A. R., Bock, J. J., Bond, J. R., Brevik, 
J. A., Contaldi, C. R., Daub, M. D., Dempsey, J. T., Goldstein, 
J. H., Holzapfel, W. L., Kuo, C. L., Lange, A. E., Lueker, M., 
Newcomb, M., Peterson, J. B., Ruhl, J., Runyan, M. C, & 
Staniszewski, Z. 2009, ApJ, 694, 1200 

Ruhl, J., Ade, P. A. R., Carlstrom, J. E., Cho, H.-M., Crawford, 
T., Dobbs, M., Greer, C. H., Halverson, N. W., Holzapfel, W. L., 
Lanting, T. M., Lee, A. T., Lerlinitch, E. M., Leong, J., Lu, W., 
Lueker, M., Mehl, J., Meyer, S. S., Mohr, J. J., Padin, S., Plagge, 
T., Pryke, C, Runyan, M. C, Schwan, D., Sharp, M. K., Spieler, 
H., Staniszewski, Z., & Stark, A. A. 2004, in Proc. SPIE, Vol. 
5498, Millimeter and Submillimeter Detectors for Astronomy II, 
ed. J. Zmuidzinas, W. S. Holland, &; S. Withington (Bellingham: 
SPIE Optical Engineering Press), 11-29 

Sayers, J., Golwala, S. R., Rossinot, P., Ade, P. A. R., Aguirre, 
J. E., Bock, J. J., Edgington, S. F., Glenn, J., Goldin, A., Haig, 

D. , Lange, A. E., Laurent, G. T., Mauskopf, P. D., & Nguyen, 

H. T. 2009, ApJ, 690, 1597 

Sehgal, N., Bode, P., Das, S., Hernandez-Monteagudo, C, 
Huffenberger, K., Lin, Y., Ostriker, J. P., & Trac, H. 2010, ApJ, 
709, 920 

Staniszewski, Z., Ade, P. A. R., Aird, K. A., Benson, B. A., Bleem, 
L. E., Carlstrom, J. E., Chang, C. L., Cho, H.-M., Crawford, 
T. M., Crites, A. T., de Haan, T., Dobbs, M. A., Halverson, 
N. W., Holder, G. P., Holzapfel, W. L., Hrubes, J. D., Joy, M., 
Keisler, R., Lanting, T. M., Lee, A. T., Leitch, E. M., Loehr, A., 
Lueker, M., McMahon, J. J., Mehl, J., Meyer, S. S., Mohr, J. J., 
Montroy, T. E., Ngeow, C.-C, Padin, S., Plagge, T., Pryke, C, 
Reichardt, C. L., Ruhl, J. E., Schaffer, K. K., Shaw, L., Shirokoff, 

E. , Spieler, H. G., Stalder, B., Stark, A. A., Vanderlinde, K., 
Vieira, J. D., Zahn, O., & Zenteno, A. 2009, ApJ, 701, 32 

Sunyaev, R. A. &; Zel'dovich, Y. B. 1972, Comments on 

Astrophysics and Space Physics, 4, 173 
Vikhlinin, A., Kravtsov, A., Forman, W., Jones, C, Markevitch, 

M., Murray, S. S., & Van Speybroeck, L. 2006, ApJ, 640, 691 
Vikhlinin, A., Kravtsov, A. V., Burenin, R. A., Ebeling, H., 

Forman, W. R., Hornstrup, A., Jones, C, Murray, S. S., Nagai, 

D., Quintana, H., & Voevodkin, A. 2009, ApJ, 692, 1060 
Vikhlinin, A., Markevitch, M., Murray, S. S., Jones, C, Forman, 

W., & Van Speybroeck, L. 2005, ApJ, 628, 655 
White, M. & Majumdar, S. 2004, ApJ, 602, 565 
Wilson, G. W., Hughes, D. H., Aretxaga, I., Ezawa, H., 

Austermann, J. E., Doyle, S., Ferrusca, D., Hernandez-Curiel, 

I. , Kawabe, R., Kitayama, T., Kohno, K., Kuboi, A., Matsuo, 
H., Mauskopf, P. D., Murakoshi, Y., Montana, A., Natarajan, 
P., Oshima, T., Ota, N., Perera, T. A., Rand, J., Scott, K. S., 
Tanaka, K., Tsuboi, M., Williams, C. C, Yamaguchi, N., & Yun, 
M. S. 2008, MNRAS, 390, 1061 

Zhang, Y.-Y., Bohringer, H., Finoguenov, A., Ikebe, Y., 
Matsushita, K., Schuecker, P., Guzzo, L., & Collins, C. A. 2006, 
A&A, 456, 55 

Zhang, Y.-Y., Finoguenov, A., Bohringer, H., Kneib, J., Smith, 
G., Kneissl, R., Okabe, N., & Dahle, H. 2008, Astronomy and 
Astrophysics, 482, 451 



16 



APPENDIX 
CLUSTER MAPS AND PROFILES 

Shown below are the SPT maps and profiles for each of the clusters in this sample. Four maps are shown for each 
cluster: the 150 GHz map (upper left), the 220 GHz map (upper right), the band subtracted map (lower left), and 
a jackknife band subtracted map (lower right). The pixel size is 0'.25, and each map is smoothed with a Gaussian 
with an FWHM of V. The 220 GHz maps contain more noise than the 150 GHz maps due to lower sensitivity and 
higher-amplitude atmospheric fluctuations in the 220 GHz band. The band subtracted map is calculated by applying a 
matched filter to the 220 GHz map and subtracting the result from the 150 GHz map. The jackknife band subtracted 
map is calculated by multiplying half of the single-observation maps by —1 before co- adding. The locations of bright 
point sources that were removed are marked by crosses. The visible structure to the east and west of the brighter 
clusters is due to the timestream polynomial subtraction (see the text), and is accounted for in the radial profiles. 

The profiles are computed from the band subtracted maps using our knowledge of the SPT beams and time-domain 
filtering. Each profile point represents the average SZ temperature decrement within a given radial bin. The dashed 
vertical line is an estimate of rsoo • The algorithm used to derive the profiles is described in Section 14.11 




Fig. Al. — A 2744 maps (left) and profile (right). The four maps are 150 GHz (upper left), 220 GHz (upper right), band subtracted 
(lower left), and jackknife (lower right). Units are /x-Kcmb- 




Fig. A2. RXCJ0217.2-5244 maps (left) and profile (right). Units are /iKcmb- 
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Fig. A3. RXCJ0232. 2-4420 maps (left) and profile (right). Units are /iK C MB- 




Fig. A4. RXCJ0336. 3-4037 maps (left) and profile (right). Units are /ii^cMB- 




Fig. A5. 



AS 0520 maps (left) and profile (right). Units are iiKqmb- 
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Fig. A6. RXCJ0528. 9-3927 maps (left) and profile (right). Units are /iK C mb- 
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Fig. A7. RXCJ0532. 9-3701 maps (left) and profile (right). Units are /ii^cMB- 
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Fig. A8. MACSJ0553. 4-3342 maps (left) and profile (right). Units are ^ C mb- 



19 





-400 r 



Fig. A9. — AS 0592 maps (left) and profile (right). Units are ^cmb- 
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Fig. A10. — A 3404 maps (left) and profile (right). Units are /ii^cMB- 
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Fig. All. — 1ES 0657-56 maps (left) and profile (right). Units are ^cmb- 
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Fig. A13. — A 3856 maps (left) and profile (right). Units are /4-Kcmb- 




Fig. A14. 



A 3888 maps (left) and profile (right). Units are /ii^CMB- 
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Fig. A15. — AS 1063 maps (left) and profile (right). Units are i^Kqmb- 



